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AB£STKACT 


A  Naval  Postgraduate  School  spacecraft  design  class 
proposed  a  multiple  beam,  phased  array,  microstrip  antenna  as 
part  of  the  preliminary  design  of  a  low  earth  orbit  communica¬ 
tion  satellite.  The  antenna  must  provide  coverage  over  the 
satellite's  entire  field  of  view  while  both  uplink  and  down¬ 
link  operate  simultaneously  on  the  same  L-band  frequency. 

This  thesis  assesses  the  feasibility  of  the  antenna 
proposed  in  that  preliminary  design.  Design  tradeoffs  for  a 
microstrip  array  constrained  by  both  available  siirface  area 
and  a  limited  mass  budget  are  examined.  Two  different 
substrate  materials  are  considered  in  terms  of  weight  and 
performance.  Microstrip  patch  theory  is  applied  to  array 
element  design  and  layout  and  antenna  array  theory  is  applied 
to  determine  phase  and  anplitude  coefficients.  The  focus  of 
the  design  is  on  obtaining  the  desired  beam  shape  and  orienta¬ 
tion,  given  euitenna  size  constraints.  A  corporate  feed  method 
is  discussed  and  a  general  design  presented. 

Antenna  performance  is  predicted  through  the  use  of  a 
computer  model  based  on  Modal  Expansion  theory  and  results  are 
plotted  in  a  series  of  graphs  which  demonstrate  the  limita- 
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I.  ZNISODOCnON 


A.  BAIC3CGR0GND 

In  J\me  1991,  Loral  Cellular  Systems  Corporation  of  Palo 
Alto,  California  submitted  a  license  request  to  the  Federal 
Coirmunications  Commission  for  authority  to  build  a  world-wide, 
itrabile  communications  system  [Ref.  1]  .  Their  proposed  system, 
called  GLOBALSTAR,  is  a  satellite  system  designed  to  provide 
Radio-Determination  Satellite  Services  (RDSS)  for  real  time 
position  location  and  tracking,  as  well  as  voice  and  data 
services.  The  system  is  intended  to  be  integrated  into  the 
existing  Public  Switched  Telephone  Network  (PSTN)  and  private 
cellular  networks  in  order  to  provide  telephone  communications 
to  users  euiywhere  in  the  world. 

In  a  project  sponsored  by  the  NASA/ttaiversities  Space 
Research  Association  Advanced  Design  Program,  students  at  the 
Naval  Postgraduate  School,  Monterey,  CA  performed  a  prelimin¬ 
ary  design  [Ref.  2]  for  a  satellite  system  based  on  the 
proposed  GLOBALSTAR  system.  The  focus  of  this  thesis  is  on 
the  design  and  performance  analysis  of  a  phased  array  L-band 
antenna,  composed  of  microstrip  patches,  which  was  proposed  as 
part  of  this  preliminary  design.  The  scope  is  limited  to  the 
design  of  the  antenna  panels  euid  microstrip  patches  and  to  the 
performance  analysis  of  overall  antenna.  While  the  beam 
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forming  network  will  be  mentioned,  a  detailed  analysis  of  the 
beam  forming  network  is  outside  the  scope  of  this  thesis. 
Similarly,  certain  mechanical  details  such  as  panel  hinging 
and  antenna  deployment  mechanisms  as  well  as  specific 
fabrication  methods  will  not  be  considered  in  detail. 

B.  SYSTEM  acMcaayr 

Using  this  system,  a  mobile  telephone  user  will  have 
access  to  world  wide  coverage  by  accessing  the  PSTN  via 
satellite.  Using  a  special  mobile  telephone  unit,  a  user 
communicates  directly  with  a  satellite  via  an  L-band  radio 
link.  The  satellite  receives  signals  from  mobile  users  and 
retransmits  these  signals  to  an  earth  based  gateway  via  an  S- 
band  trunk  link.  In  the  reverse  process,  signals  from  the 
gateways  are  received  by  the  satellite  via  the  S-band  link  and 
then  retransmitted  to  mobile  users  via  the  L-band  link. 
Gateways  are  distributed  across  the  globe  and  are  connected  to 
the  PSTN. 

Acquisition,  synchronization,  and  satellite  beam  hand  off 
coordination  is  provided  for  the  gateways  by  a  Satellite 
Operation  Control  Center  (SOCC)  as  depicted  in  Figure  1.1. 
The  information  provided  by  the  SOCC  allows  each  gateway 
station  to  track  and  communicate  with  the  satellites  in  its 
field  of  view.  Ephemeris  data  is  provided  to  the  SOCC  by 
several  Tracking  Telemetry  and  Control  Stations  (Tr&C)  which 
obtain  the  data  through  satellite  obseirvations .  Overall 
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communications  network  management  is  provided  by  a  Network 
Control  Center  (NCC) . 


OATEWMirS 


Figure  1.1  System  Architecture 


C.  SYSTQl  DESIGN 

An  orbital  altitude  of  750  nm  gives  a  satellite  field  of 
view  of  approximately  108  degrees  (assianing  a  minimum 
elevation  angle  of  10  degrees  for  the  satellite  observer) . 
This  allows  a  constellation  of  48  identical  satellites  in 
circular  orbits  to  provide  continuous  coverage  between  75 
degrees  North  latitude  and  75  degrees  South  latitude. 
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The  constellation  is  divided  into  eight  orbital  planes  at 
an  interval  of  30  degrees  with  six  satellites  in  each  plane. 
There  is  no  provision  for  on-orbit  spares.  Communication 
between  the  satellite  and  gateways  is  accomplished  via  an 
uplink  and  downlink  at  center  frequencies  of  6533.25  MHz  euid 
5207.75  MHz  respectively,  each  with  a  bandwidth  of  16.5  MHz. 
A  mobile  user  communicates  with  a  satellite  via  an  L-band  link 
operating  at  a  center  frequency  of  1618.25  MHz  and  a  bandwidth 
of  16.5  MHz.  A  summary  of  the  link  budget  as  specified  by 
[Ref.  1]  is  contained  in  Appendix  A. 

Each  satellite  provides  ccxnmunications  coverage  over  its 
entire  field  of  view  through  the  use  of  six  non-overlapping 
elliptically  shaped  spot  beams.  To  maximize  the  amoimt  of 
time  a  mobile  user  is  illuminated  by  any  one  particular  beam, 
beams  are  oriented  such  that  their  major  axes  are  parallel  to 
the  spacecraft  velocity  vector.  Fig\ire  1.2  shows  a  typical 
coverage  pattern  for  one  satellite. 

In  order  to  minimize  interference  between  channels,  the 
six  beams  operate  in  pairs  such  that  maximum  beam  separation 
is  achieved  on  the  ground.  This  is  accomplished  by  switching 
beams  on  and  off  in  pairs  so  that  opposite  pairs  are  active 
for  10  ms  out  of  a  60  ms  duty  cycle.  During  the  duty  cycle, 
beam  pairs  alternate  between  transmit  and  receive  functions. 
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D.  ANTQ1N21  PERFORMANCE  REQDIREMEMTS 


The  L-band  antenna  for  this  system  must  not  only  provide 
adequate  gain,  but  must  also  correctly  shape  each  of  the  six 
spot  beams  while  maintaining  low  sidelobes.  The  antenna  must 
also  allow  the  communications  payload  to  maintain  a  constant 
bit  error  rate  (BER)  .  Bit  error  rate  is  a  fmction  of  carrier 
to  noise  rauio  (C/N)  and  C/N  is  a  function  of  path  length. 
Because  antenna  beams  must  be  statically  scanned  off  center  to 
cover  the  FOV,  the  path  length  for  each  of  the  six  antenna 
beams  is  different.  Path  lengths  vary  from  a  minimum  of  1748 
km  (measured  along  antenna  boresight)  for  the  inner  most 
beams,  to  a  maximum  of  9066  km  for  the  outer  most  beams.  In 
order  to  maintain  a  constant  BER,  an  isofliix  antenna  is 
required.  Table  I-l  summarizes  the  L-band  antenna  performance 
parameters . 


TABLE  I-l.  hNTBSftOi  PARAMETERS 


f  0=1618. 25 
BW=16 . 5 

(MHz) 

Beam 

1  &  6 

Beam 

2  &  5 

Beam 

3  &  4 

Center 

45* 

30* 

10* 

BW  (-4  dB) 

25*  X  110* 

35*  X  120* 

40*  X  120* 

BW  (total) 

40*  X  120* 

50*  X  120* 

70*  X  140* 

Gain  (dBi) 

6.5 

5.5 

1 

o 

• 

Note:  Satellite  FOV  is  108* 


II.  DESIGN  C9GNSIDERATI0NS 


A.  SPACECRAFT  OONETGDRATICN 

1.  Overview 

'The  spacecraft  is  designed  to  be  placed  into  orbit 
using  the  Delta  7925  launch  vehicle.  A  Satellite  Launch 
Dispenser  (SLD)  was  designed  so  that  a  single  laimch  vehicle 
can  place  an  entire  plane  of  satellites  into  orbit  at  one 
tinte.  The  SLD  is  essentially  a  chest  of  drawers  where  the 
spacecraft  are  the  drawers.  The  spacecraft  slide  into  the  SLD 
along  rails  and  conpress  a  spring  when  properly  seated.  Bach 
spacecraft  is  then  held  in  place  by  explosive  bolts.  At  the 
appropriate  time  during  orbit  insertion,  the  explosive  bolts 
are  fired  and  a  spacecraft  is  ejected  into  the  proper  orbit. 

2.  Size  CoKistraints 

The  Delta  payload  shroud  and  the  SLD  impose  serious 
size  constraints  on  overall  spacecraft  size.  In  order  to  fit 
into  the  space  provided  by  the  SLD,  the  spacecraft  dimensions 
can  not  exceed  1.4  x  1.4  x  0.6  meters.  The  space  allowed  for 
the  L-band  antenna  alone  is  1.4  x  1.4  x  0.1  meters  in  the 
stowed  configuration.  The  preliminary  design  allowed  an 
antenna  mass  budget  of  10  kg  however  this  proves  to  be 
unrealistic . 


7 


B.  ANTENNA  CJONFIGORATION 

In  order  to  achieve  beam  steering  and  beam  shaping 
characteristics  necessary  to  meet  the  antenna  performance 
requirements,  a  planar  phased  array  antenna  similar  to  the  L- 
band  antenna  of  [Ref.  3]  was  considered.  A  planar  arxay 
allows  both  beam  shaping  and  beam  steering  in  any  direction. 
For  sinplicity,  rectangular  microstrip  patch  array  elements 
were  chosen.  A  microstrip  antenna  is  both  lightweight  and 
easily  fabricated  using  common  printed  circuit  board  etching 
techniques  [Ref.  4] .  Disadvantages  usually  associated  with 
microstrip  antennas,  specifically  large  surface  area  require¬ 
ments  and  narrow  bandwidth  [Ref.  5]  ,  were  foimd  to  be  insigni¬ 
ficant  in  this  application.  Folding  array  panels  are  employed 
so  that  the  antenna  can  be  stowed  in  the  space  available  and 
conveniently  deployed  via  a  spring-loaded  mechanism  once  the 
spacecraft  is  in  the  proper  orbit.  Each  array  panel 
corresponds  to  one  of  the  six  antennas  necessary  to  form  the 
appropriate  beam.  Figure  2.1  shows  the  phased  array  antenna 
in  the  deployed  configuration. 

C.  MZC310SIRIP  PATCH  THEORY 

The  simplest  and  most  common  type  of  microstrip  antenna  is 
the  rectangular  patch  mounted  on  a  dielectric  substrate  of 
thickness  h,  shown  in  Figure  2.2. 


8 


Figure  2.2.  Rectaxigular  Microstrip  Patch 


1.  Radiatloa  Mechanisms 

The  radiation  mechanism  of  a  microstrip  patch  can  be 
modeled  in  one  of  several  ways,  each  with  its  advantages  and 
disadvantages . 

a.  The  Wire  Grid  Model 

A  microstrip  patch  can  be  modeled  as  a  grid  of 
fine  wires  [Ref.  6]  .  The  cnirrents  on  the  wires  are  solved  for 
numerically  and  when  found,  provide  accxirate  information  about 
the  antenna  pattern  and  input  inpedcuice.  The  disadvantage  of 
this  method  is  that  it  requires  an  extraordinary  amount  of 
computer  resource  [Ref.  7]  . 
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b.  The  Modal  E}9>azisiGa  Model 


The  modal  expansion  model,  which  is  valid  for 
electrically  thin  (h«l)  substrates,  represents  the  antenna  as 
a  thin  TM  mode  cavity  boimded  by  magnetic  walls  [Ref.  5] 
within  which  the  dominant  TM  mode  is  excited  along  with  other 
non- resonant  modes.  As  long  as  the  dominant  mode  is  suffi¬ 
ciently  excited,  this  model  provides  accurate  performance 
information  [Ref.  7],  including  the  effects  of  feed  location. 
The  ability  to  accoxmt  for  a  feed  probe  and  handle  complex 
patch  shapes  are  the  chief  advantages  of  this  model, 
c.  Hie  Transmission  Line  Model 

The  transmission  line  model  considers  the  patch  to 
be  a  line  resonator  with  radiation  due  to  fringing  fields  at 
the  open  circuited  ends  [Ref.  6]  .  The  radiator  consists  of 
two  radiating  slots  separated  by  a  distance  L.  The  slots  are 
formed  at  the  edges  of  the  patch  element,  between  the  patch 
and  the  ground  plcuie  below.  Field  variations  along  the 
radiating  edges  are  ignored.  The  chief  disadvantage  of  the 
treuismission  line  model  is  that  it  is  only  valid  for  rectan¬ 
gular  patches.  The  input  impedance  of  the  transmission  line 
model  is  also  highly  dependent  on  feed  location  and  fails  to 
accxirately  predict  input  impedance  for  all  feed  locations 
[Ref.  7]  .  However,  formulas  have  been  developed  to  predict 
input  impedance  for  any  feed  location  [Ref.  6]  . 

Despite  its  dratdjacks,  the  transmission  line 
model's  main  advantage  is,  that  for  rectangular  patches,  it 
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provides  single  design  formulas  which  give  reasonably  accurate 
performance  predictions. 

2 .  Patch  Parameters 

The  dimensions  of  a  rectangular  patch,  as  well  as 
bandwidth  and  gain,  are  determined  by  the  operating  frequency 
of  the  antenna,  the  relative  dielectric  constant,  and 
thickness  of  the  substrate  material.  The  following  formulas 
are  based  on  the  transmission  line  model, 
a.  Width  and  Lesufth 

The  width  and  length  of  a  rectangular  microstrip 
patch  are  given  by: 


(cm) 


(2.1) 


L  =  2 —  -  2Ai  (cm) 


where 

c  =  speed  of  light  (m/s) , 

=  operating  frequency  (MHz) , 

=  relative  dielectric  constant, 

=  effective  dielectric  constant. 


e 


e^+1 ^  e^-1 
2  2 


12h\-i 

W  ) 


(2.2) 


(2.3) 
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(2.4) 


(eg+0.3)  (-^+0.264) 

Ai  =  0.412A  - - 

(e  -0.258)  (-^+0.8) 
h 

and 

h  =  dielectric  thickness  (cm)  [Ref.  6] . 

b.  Bandsd-dth 

The  bandwidth  of  a  microstrip  antenna  is  defined 
in  terms  of  the  antenna's  quality  factor  (Q)  as  follows: 

BW  =  (2.5) 

O^vSm 

Where  VSWR  is  less  than  a  specified  value  (2:1  or  1.5:1  for 
exanple)  and  VSWR=1  at  the  operating  frequency  [Refs.  5  and 
6]  ,  the  bandwidth  is  dependent  on  both  the  relative  dielectric 
constant,  e^,  and  thickness  of  the  substrate.  Thicker  sub¬ 
strates  and  lower  values  of  give  larger  bandwidths. 
Typically,  microstrip  antennas  have  bandwidths  on  the  order  of 
a  few  percent  of  the  operating  frequency. 

c.  Gain 

Gain  is  also  affected  by  svibstrate  thickness  and 
relative  dielectric  constant.  Gain  is  inversely  proportional 
to  and  directly  proportional  to  substrate  thickness.  To 
achieve  the  necessary  gain,  microstrip  patches  can  be  cunranged 
to  form  lineou:  and  planar  arrays. 
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D.  ARRAYS 


The  radiation  pattern  of  a  single  array  element  is 
relatively  wide  and  typically  has  low  gain.  In  order  to 
achieve  higher  directivity  and  additional  gain,  antenna 
elements  can  be  arranged  to  form  linear  or  planar  arrays. 

1.  Theory 

The  total  field  of  an  array  is  given  by  the  vector 
addition  of  the  fields  radiated  by  individual  elements  [Ref. 
8]  .  The  required  directivity  is  achieved  by  arranging 
individual  element  excitations  such  that  fields  produced  by 
each  element  interfere  constructively  in  the  desired  direction 
and  interfere  destructively  in  all  other  directions.  The 
radiation  pattern  of  an  array  is  obtained  by  multiplying  the 
field  produced  by  a  single  element  by  an  array  factor  (AF) . 
The  array  factor  for  an  array  of  N  identical  elements  is  given 
by  [Ref.  8]  as 


w 


AF  =  5^ 


(2.6) 


where 

N  =  nxamber  of  elements, 

I„  =  element  excitation  coefficient. 
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ijr  =  Jed  X  cos  (0)  +P  , 


(2.7) 


and 

k  =  wave  number, 

d  =  distance  between  elements  (m) , 

0  =  direction  of  main  beam  (degrees) , 
p  =  phase  shift  between  elements  (degrees)  [Ref.  8] . 


The  main  beam  can  be  scanned  in  any  direction  by 
adjusting  the  progressive  phase  shift  between  elements. 

M  arrays  of  N  elements  can  be  arranged  to  form  a 
planar  array  which  will  allow  the  beam  to  be  scanned  in  any 
direction.  Additionally,  array  elements  can  be  excited 
non-uniformly  to  control  the  beam  shape.  If  a  coordinate 
system  such  as  that  of  Figure  2.3  is  chosen,  it  is  shown  in 
[Ref.  8]  that  the  array  factor  in  the  X  direction  is  given  by 


(Jedj^sin (8)  COB  (4)  ♦Px) 


(2.8) 


iii«i 


Similarly,  the  array  factor  in  the  Y  direction  is  given  by 

N 


j(n-l)  (lccr^in(e)sin(^)  ♦Py) 


(2.9) 


tt*l 


where 


=  excitation  coefficient  for  m  elements  and 
J,  =  excitation  coefficient  for  n  elements. 
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Figure  2.3.  Array  Coordisate  System 

The  total  array  factor  is  the  product  of  and  AF^ 
and  is  given  by 


M  n 

I  gj(w-l)  (Wj^in(e)c08(^).pj^  J  (lc</y8in(a)8in(4)  ♦Pj,) 

in»l 


N 

c 

n«l 


(2.10) 

2 .  Beam  Synthesis 

The  becun  in  a  phased  array  antenna  is  shaped  by  non- 
uniform  excitation  of  the  array  elements .  This  non-uniformity 
Ceui  be  in  amplitude,  phase,  or  both.  Using  general  synthesis 
[Ref.  9]  and  optimization  [Refs.  10  and  11]  techniques,  beams 
of  almost  any  shape  and  contour  can  be  formed.  For  less 
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stringent  beam  shape  criteria,  a  simpler  synthesis  method 
based  on  the  Foxirier  transform  method  presented  in  [Ref.  8] 
has  been  derived. 

In  the  Fourier  transform  method,  excitation 
coefficients  are  related  to  the  Fourier  transform  of  the  array 
factor  expressed  as  a  function  of  \|r,  which  is  itself  a 
function  of  0  and  4>.  Excitation  coefficients  are  given  by 


=  -^/_*AF(ilr)e^«+  dijr  , 
-N  n  ^  N 


(2.11) 


for  an  odd  niimber  of  elements,  and 


1  -if  1* 


(2.12) 


-M  ^  m  ^  -1 


and 


In,  = 


-j[ 


di|r 


1  i  m  i  Af 


(2.13) 


for  an  even  number  of  elements  [Ref.  8] . 

Equations  2.11  -  2.13  yield  conplex  numbers  which 
represent  the  amplitude  and  phase  excitation  of  an  element. 
Ihe  total  excitation  coefficient  is  obtained  as  the  product  of 
and  I.. 

n  m 


17 


III.  ANIEtOOl  DESIGN 


A.  MICROSTSIP  PATCH  DESIGN 
1.  Substarate  Material 

The  first  step  in  the  design  process  is  the  selection 
of  a  dielectric  substrate  material.  The  dominant  feattires  of 
a  microstrip  array  are  controlled  by  sxibstrate  parameters  such 
as  thickness  and  permittivity  more  than  by  the  particular 
element  type  [Ref.  12] . 

Because  weight  is  a  primary  consideration  in 
spacecraft  design,  a  substrate  must  be  chosen  which  not  only 
has  satisfactory  dielectric  properties,  but  which  also  has  low 
density.  There  are  a  variety  of  substrate  materials  available 
with  a  wide  range  of  relative  dielectric  constants  and 
densities.  Tables  sumnarizing  various  materials  and  their 
relevant  properties  can  be  found  in  [Refs.  5  and  6]  . 

Two  representative  substrate  materials  which  were 
considered  for  this  design  are  Duroid  (RT  5870) ,  produced  by 
Rogers  Corporation  of  Cheuidler,  AZ,  and  Kevlar  (HI  4093) , 
produced  by  Arlon  Electronic  Substrates  Division  of  Rancho 
Cucamonga,  CA. 

Studies  have  shown  [Ref.  13]  that  efficiency  and 
bandwidth  are  nearly  independent  of  patch  shape  (i.e., 
circular  vs  rectangular)  and  are  determined  mainly  by 
substrate  thickness  and  permittivity. 


In  order  to  achieve  good  resolution  in  beam  shape  and 
steering,  the  number  of  array  elements  must  be  maximized 
within  available  array  area  and  element  spacing  constraints. 
This  is  acconplished  by  madcing  the  microstrip  patches  as  small 
as  possible  while  maintaining  acceptable  performance  charac¬ 
teristics.  For  a  given  operating  frequency,  patch  size  is 
determined  by  relative  dielectric  constant  and  stabstrate 
thickness .  Figure  3 . 1  shows  the  effect  of  relative  dielectric 
on  patch  width  based  on  equation  (2.1)  .  Figure  3.2  shows  the 
effect  of  dielectric  thickness  (h)  on  patch  length  for  a  given 
operating  frequency  and  relative  dielectric  constant,  based  on 
equation  (2.2). 

Duroid  substrates  are  available  with  ranging  from 
2.1  to  10.0.  Higher  yields  a  smaller  patch  size  at  the 
expense  of  bandwidth  and  gain.  Rogers  Duroid  RT  5870  has 
equal  to  2.23  which  allows  higher  efficiency  and  bandwidth 
performance  but  the  density  of  RT  5870  Duroid  is  2.2  gm/cm?. 
Arlon's  HI  4093  Kevlar  material  on  the  other  hand,  has  a 
relative  dielectric  constant  of  3.9  and  a  density  of  1.38 
gm/cm?. 

For  a  given  operating  frequency  and  dielectric 
thickness,  Kevlar  yields  a  smaller  patch  size  than  Duroid  and 
although  efficiency  and  bandwidth  performance  are  reduced,  the 
results  are  still  within  acceptable  parameters.  Therefore, 
the  lower  mass  of  Kevlar  far  outweighs  its  disadvantages. 
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Figure  3.1.  Effect  of  on  Width  of  a  Microstrip  Patch 


Figure  3.2 
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Using  a  solid  Duroid  dielectric  0.3  cm  thick,  the  mass 
of  the  Slabs trate  alone  would  be  36.7  kg  (for  an  antenna  area 
ofl.4mx4.0m).A  Kevlar  substrate  with  the  same  dimensions 
would  yield  an  antenna  mass  of  23.3  kg.  Appendix  B  contains 
detailed  mass  estimates  for  antenna  components. 

2.  Patch  DimeiDsioins 

a.  Thickness 

Because  the  antenna  substrate  and  ground  plane 
must  form  the  array  panels,  they  must  be  made  thick  enough  to 
provide  the  required  rigidity.  The  array  panels  must  be  rigid 
enough  to  withstand  structural  loads  of  launch  and  deployment 
as  well  as  resist  warpage  due  to  thermal  cycling.  A  siabstrate 
thickness  of  at  least  0.3  cm  was  selected  and  assumed  adequate 
to  meet  these  criteria.  The  antenna  mass  could  be  reduced 
without  significantly  affecting  electrical  and  mechanical 
performance  by  supporting  a  thin  dielectric  substrate  on  a 
honey-comb  structure  such  as  that  described  in  [Refs.  3  and 
4] .  The  effect  of  raising  the  height  of  the  patch  above  the 
ground  plane  is  to  increase  bandwidth  [Ref.  13],  while 
reducing  dielectric  thickness  would  have  the  opposite  effect . 

b.  Area 

The  equations  presented  in  Cliapter  II  can  be  used 
as  a  starting  point  for  determining  patch  dimensions  which  can 
later  be  optimized  using  a  more  sophisticated  model.  For  an 
operating  frequency  of  f  ^=1618. 25  MHz  and  a  relative 
dielectric  constant  e^=3.9,  equations  (2.1) -(2.4)  yield  a 
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patch  length  of  5.8  cm  and  a  width  of  7.2  cm  (note  that  the 
width  is  defined  as  the  long  dimension  of  the  rectangle)  .  The 
gain  of  a  single  patch  is  approximately  5  dB.  Detailed  cal¬ 
culations  for  length,  width,  efficiency,  directivity,  and  gain 
calculations  for  a  single  patch  are  contained  in  i^pendix  C. 

A  FORTRAN  computer  program  which  makes  use  of 
modal  e3q>ansion  theory  was  used  in  an  iterative  process  to 
optimize  the  values  for  patch  length  given  a  constant  width. 
A  listing  of  the  program  is  contained  in  Appendix  F.  The 
final  values  for  length  and  width  are  summarized  in  Table  III- 
1  along  with  the  operating  parameters  for  a  single  microstrip 
patch. 


TABLE  III-l.  PARAMETERS  FOR  A  SINGLE  MICROSTRIP  PATCH 


e  *3.9,  h=0.3 
(cm) 

f, =1618. 25  MHz 

L  (cm) 

W  (cm) 

6  (dB) 

%  BW 

Transmission  Line 

4.6 

5.9 

5.02 

2.42 

Modal  Expansion 

4.44 

5.9 

5.02 

2.62 

B.  PLANAR  ARRAY  DESIGN 

Six  array  panels  are  needed  to  produce  the  six  ellipti- 
cally  shaped  beams.  The  beams  are  symmetrical  from  outboard 
to  inboard  with  beams  1  &  6  being  identical,  beams  2  &  5 
identical,  and  beams  3  &  4  identical.  The  total  area  avail¬ 
able  for  each  array  panel  in  the  deployed  configuration  is 
0.9338  m^. 
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1.  Element  Spacing 

As  previously  stated,  the  number  of  array  elements 
needs  to  be  maximized  in  order  to  achieve  good  beam 
resolution.  At  the  same  time,  the  array  elements  must  be 
placed  far  enough  apart  to  avoid  mutual  coupling  which  would 
degrade  performance  by  increasing  side  lobe  levels  (SLL)  .  It 
has  been  shown  [Ref.  14]  that  for  an  element  spacing  of  d  > 
X/2,  mutual  coupling  will  not  cause  significant  side  lobe 
levels.  Measurements  made  on  L-band  microstrip  arrays  have 
shown  that  mutual  coupling  levels  of  -25  dB  or  less  can  be 
achieved  with  proper  element  spacing  [Refs.  14  and  15]  .  If 
the  elements  are  separated  by  much  more  than  one  half 
wavelength,  imwanted  grating  lobes  are  produced.  An  element 
spacing  of  about  0.55X,  has  been  shown  to  give  good  results 
[Ref.  3]  .  To  avoid  grating  lobes  and  mutual  coupling  effects, 
the  patch  spacing  for  this  design  was  chosen  as  0.1017  m 
(0.55X,).  This  spacing  allows  a  4  x  9  element  array  to  be 
placed  on  each  of  the  six  antenna  panels. 

2 .  Excltatioga  Ooef  f  iclents 

The  Fourier  treuisform  method  which  was  previously 
developed  for  a  planar  array  was  used  to  obtain  element  phase 
and  anplitude  excitation  coefficients.  Array  coefficients  for 
a  single  beam  are  shown  in  Table  II 1-2.  The  resulting  array 
factors  are  plotted  in  Figure  3.3  and  Figure  3.4.  The  excita¬ 
tion  coefficients  are  tabulated  in  J^>pendix  D  and  array  factor 
plots  for  each  beam  are  presented  in  7^>pendix  E. 
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TABLE  IZI-2.  BEAM  1  EXdTATIGN  COEFFICIENTS 


Element 

Nimber 

1 

2 

3 

4 

1 

0.013 
^52. 36* 

0.046 

4.-127.6* 

0.046 

4.-127.6* 

0.013 

4.52.36* 

2 

0.019 

^-57.73* 

0.066 

4.129.3* 

0.066 

4.129.3* 

0.019 

4.-57.73* 

3 

0.024 

A.-153.8‘ 

0.083 

/IL28.18* 

0.083 

4.28.18* 

0.024 

4.-153.8* 

4 

0.027 
^103 . 1* 

0.094 
4.-76. 9* 

0.094 
4.-76. 9* 

0.027 

4.103.1* 

5 

0.028 

AO.O* 

0.098 

4.180.0* 

0.098 

4.-180.0* 

0.028 

4.0.0* 

6 

0.027 

A-103.1’ 

0.094 

4.76 . 9* 

0.094 

4.76 . 9* 

0.027 

4.-103.1* 

7 

0.024 
L152 . 8* 

0.083 

4.-28.18* 

0.083 

4.-28.18* 

0.024 
4.153 . 8* 

8 

0.019 
LSI. 12' 

0.066 

4.-129.3* 

0.066 

4.-129.3* 

0.019 

4.57.73* 

9 

0.013 

^-52.36* 

0.046 

4l127.6* 

0.046 

4.127.6* 

0.013 

4.-52.36* 

3 .  Antenna  Gain 

The  inaxinTum  array  gain  is  the  s\jm  of  individual 
element  gains  less  losses  associated  with  the  array  feed.  The 
maximum  gain  is  20.6  dBi.  This  is  well  above  the  maximum 
required  (6  dBi)  for  any  euitenna  panel,  however,  a  significant 
amount  of  loss  will  occur  in  the  feed  structure  and  must  be 
offset  by  the  excess  gain. 


Figure  3.3.  Beam  1  Array  Factor  (4>b90*) 
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4.  Feed  Methods 


There  are  a  variety  of  feed  methods  for  microstrip 
patch  arrays  and  for  each  method  there  exists  many  variations . 
In  general,  feed  methods  are  grouped  into  three  categories  - 
probe,  corporate,  and  triplate  feeds  -  and  are  briefly 
discussed  next. 

a.  Probe  Feed 

In  the  probe  feed  method,  each  patch  is  connected 
to  the  feed  network  via  a  coaxial -fed  probe  which  extends 
through  the  ground  plane  and  substrate.  The  probe  is  physi¬ 
cally  attached  to  the  patch  by  a  solder  joint.  Impedance 
matching  between  the  coeix  and  array  element  can  be  accom¬ 
plished  for  thick  substrates  by  designing  a  teardrop- shaped 
probe  which  cancels  out  probe  inductance  [Ref.  3]  or  by 
varying  the  diameter  of  the  probe  for  thin  substrates  [Ref. 
16]  .  Obviously  it  would  be  inpractical  to  construct  an  array 
where  each  element  had  its  own  coaxial  feed  cable.  The  weight 
and  space  requirements  would  be  prohibitive  for  spacecraft 
applications . 

b.  CSorpoFate  Feed 

A  corporate  feed  offers  a  practical  solution  to 
the  array  feed  problem  by  allowing  one  coaxial  line  to  feed  a 
network  composed  of  microstrip  transmission  lines.  A  series 
of  power  dividers  delivers  power  to  and  from  the  radiating 
elements.  The  microstrip  transmission  line  provides  equal 
path  length  to  all  elements  maintaining  phase  between 
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elements.  The  two  major  disadvantages  of  the  corporate  feed 
are  that  there  is  not  enough  room  on  a  planar  arxay  substrate 
for  radiating  elements  and  the  feed  network,  and  high  line 
loss  occurs  due  to  long  path  lengths.  The  second  disadvantage 
can  be  overcome  by  integrating  active  elements  into  the 
substrate  [Refs.  17  and  18]  or  by  breaking  the  feed  network 
into  sxibsections,  each  of  which  is  fed  by  a  low  loss  coaxial 
cable  [Ref.  19] . 

c .  nriplate  Feed 

In  a  triplatcr  feed  (also  called  a  hybrid  feed) , 
the  feed  network  is  etched  on  its  own  substrate  layer  and 
either  laminated  to  the  underside  of  the  radiating  element's 
sxibstrate  or  attached  to  the  backside  of  the  ground  plane. 
The  former  method  requires  the  radiation  of  the  feed  network 
to  be  accounted  for  in  the  antenna  design  while  the  latter 
method  eliminates  feed  radiation  effects  and  divides  the 
design  problem  into  two  parts:  the  array  itself  and  the  feed 
network . 

(1)  Feed  Network  in  Front  of  Ground  Plane.  A 
triplate  feed  located  in  front  of  the  groxand  plane  and 
directly  vutder  the  radiating  element's  substrate  simplifies 
the  array  construction  at  the  expense  of  increased  design 
complexity.  The  metal  pins  used  to  attach  the  array  element 
to  the  feed  network  are  eliminated  through  electromagnetic 
coupling.  Impedance  matching  is  controlled  by  the  amoaont  of 
overlap  between  the  microstrip  patch  and  the  feed  network. 
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Radiation  effects  can  be  largely  reduced  by  locating  the  feed 
network  as  close  to  the  grovind  plane  as  possible.  [Ref.  13] 
(2)  Triplate  Behind  Ground  Plane.  If  the 
corporate  feed  network  is  located  behind  the  ground  plane, 
ra  ’ation  effects  from  the  feed  network  are  eliminated  at  the 
e:  use  of  a  more  conplicated  and  expensive  physical  construc¬ 

tion  involving  plated  through  holes  and  connecting  pins. 

5 .  Feed  Design 

Although  antenna  feed  design  is  discussed  in  [Ref. 
20],  a  complete  design  is  outside  the  scope  of  this  thesis. 
A  brief  discussion  of  a  possible  feed  method  is  included  to 
prove  the  feasibility  of  providing  the  required  array  element 
excitations . 

A  triplate  corporate  feed  located  behind  the  antenna 
ground  plane  could  be  xised.  This  feed  method,  based  on  those 
presented  in  [Refs .  6  and  19] ,  would  simplify  antenna  perform¬ 
ance  analysis  by  eliminating  feed  radiation  effects.  Since 
the  antenna  is  a  statically  scazmed  array,  the  required  phase 
shifts  and  inpedance  matching  could  be  built  into  the 
microstrip  feed  lines  using  microstrip  transmission  line 
techniques.  Figure  3.5  shows  one  possible  feed  layout. 

C.  SOMMABY 

An  exploded  view  of  one  antenna  panel  is  illustrated  in 
Figure  3.6. 
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IV.  PERFORMANCE  ANALYSIS 


A.  ANTENNA  MODEL 

Chapter  II  described  the  wire  grid,  transmission  line,  and 
modal  expansion  models  which  are  used  to  analyze  and  predict 
the  performance  of  microstrip  antennas.  The  wire  grid  model 
is  unsuitable  for  use  on  any  but  the  simplest  of  problems  due 
to  the  conputational  resoiirce  required.  Because  the  trans¬ 
mission  line  model  and  the  nKxial  expansion  model  do  not 
account  for  surface  waves  in  the  dielectric  substrate  and  do 
not  accovint  for  mutiial  coupling,  they  do  not  provide  the  most 
accurate  solution  available  [Ref.  16] .  However,  it  will  be 
shown  that  for  arrays  with  proper  element  spacing,  and  all  but 
extreme  scan  angles,  the  modal  e3q)ansion  modal  provides 
adequate  results  with  relative  ease  of  confutation.  The  areas 
where  the  modal  expansion  model  fails  will  be  discussed  and 
shown  to  be  insignificant  for  this  analysis. 

1.  Infinite  vs.  Finite  Arrays 

The  analysis  of  an  infinite  array  yields  adequate 
results  for  the  inner  array  elements  and  accoxonts  for  a  change 
in  inpedance  with  varying  scan  angle .  Infinite  array  euialysis 
ignores  edge  effects  due  to  surface  waves  in  the  dielectric. 
The  analysis  of  a  finite  array  includes  these  edge  effects  but 
at  the  esqxense  of  more  difficult  calculations.  It  has  been 
shown  that  surface  wave  power  increases  with  dielectric 
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thickness  and  that  for  a  thin  dielectric  (h  <  0.03X)  ,  surface 
wave  power  is  negligible  [Ref.  14]  .  It  has  also  been  shown 
that  surface  wave  power  diminishes  with  increasing  array  size 
and  that  for  array  sizes  of  7  x  7  or  greater  an  infinite 
approximation  is  acceptable  [Ref.  21] . 

2 .  Edge  Effects 

Diffractions  cause  by  the  finite  edge  of  the  ground 
plane  play  a  significant  part  in  the  radiation  pattern  at  wide 
scan  angles  and  in  the  backlobe  region.  The  Geometric  Theory 
of  Diffraction  (GTD)  can  be  combined  with  other  cuialysis 
methods  to  provide  an  accurate  antenna  performance  model  in 
these  regions.  Diffraction  edge  effects  can  be  ignored  for 
sccui  angles  less  than  about  70*  [Ref.  22]  . 

3 .  Jxxpat  Impedance 

Although  the  modal  expansion  model  enables  calculation 
of  input  inpedance,  patch  currents  near  the  feed  point  are  not 
modeled.  This  results  in  inaccxirate  impedance.  It  has  been 
shown  that  for  substrate  thicknesses  of  less  than  about  0.02X., 
the  idealized  feed  model  provides  adequate  results  [Ref.  23]  . 

B.  SOFTWARE  MODELS 

1.  Method  of  Moments 

The  method  of  moments  solves  for  \mknown  cxirrents  on 
the  siirface  of  the  microstrip  patch.  If  a  dielectric  is 
present,  the  dielectric  Green's  function  is  used  to  solve  the 
electric  field  integral  equation.  This  method  automatically 
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accoimts  for  mutual  coupling  and  provides  accurate  results  for 
antenna  performance.  The  method  of  moments  depends  on  the 
accurate  calculation  of  the  elements  in  the  in^dance  matrix 
for  its  precision  [Ref.  24]  eind  requires  a  large  amount  of 
cotrputer  resovurce .  Several  methods  have  been  developed  to 
calculate  method  of  moment  solutions  more  efficiently,  [Refs. 
25  and  26] ,  but  comparisons  with  results  from  other  models  do 
not  show  significant  differences  [Ref.  24] . 

2.  Modal  E}q>ansian  Model 

The  modal  expansion  model  was  discussed  in  Chapter  II . 
The  limitations  of  this  model  have  been  shown  in  previous 
discussions  to  be  insignificant  for  this  application.  Because 
of  its  simplicity,  the  modal  expcuision  model  was  used  to 
predict  the  performance  of  this  antenna.  The  software  program 
used  for  this  analysis  was  written  in  FORTRAN  and  obtained 
from  the  Jet  Propulsion  Laboratory  in  Pasadena,  CA.  A  sovirce 
code  listing  for  this  program  is  contained  in  ;^pendix  F. 

C.  MDOELIMG  SCBARIO 

Using  the  excitation  coefficients  previously  obtained, 
antenna  patterns  for  each  of  the  six  beams  were  generated  for 
scan  angles  of  0=-9O®  to  +90®,  <j>=0  and  0=-9O®  to  +90®,  <|>=90. 
The  results  were  plotted  in  order  to  asses  whether  or  not  each 
beam  met  the  design  specification  in  terms  of  scan  angle  and 
beam  width.  A  contour  plot  of  the  approximate  groxmd 
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footprint  of  each  beam  was  also  generated  in  order  to  assess 
antenna  groimd  coverage. 

D.  KESOLTS 

1.  RauUaticn  Patterns 

The  radiation  pattern  obtained  for  beam  1  is  plotted 
in  Figiire  4.1.  Beam  patterns  for  the  remaining  beams  are 
contained  in  Appendix  G. 

The  model  restilts  demonstrate  that  both  beam  shape  and 
beam  orientation  are  attainable  for  beam  1.  Similar  results 
were  obtained  for  beams  2  through  6  with  one  notable 
exception.  The  specification  called  for  beams  2  and  5  to  be 
scanned  30  degrees  from  center.  When  the  beams  were  steered 
to  exactly  -30  degrees  and  30  degrees  respectively,  the  result 
resembled  a  broad  side  array  beam  pattern.  Examination  of  the 
antenna's  normal  beam  pattern  reveals  a  null  which  prevents 
this  antenna  from  being  scanned  to  exactly  30  degrees.  By 
adjusting  the  scan  euigle  to  -27  degrees  for  beam  2  and  27 
degrees  for  beam  5,  acceptable  beam  patterns  were  obtained. 
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Figure  4.1.  Beam  1  Radiation  Pattern  (^at90') 


2 .  Pattern  Separation 

Recall  that  the  beams  cure  required  to  operate  in 
transmit /receive  pairs  such  that  the  satellite  uplink  and 
downlink  can  be  operated  simultaneously  at  the  same  £rec[uency. 
This  requires  adequate  pattern  separation  to  avoid  mutvial 
interference  at  the  eeu:th's  surface.  Figture  4.2  demonstrates 
that  this  design  achieves  adequate  pattern  separation  for 
beams  1  and  4.  Similar  resxilts  were  obtained  for  beam  pair  2 
and  5  and  pair  3  and  6  and  are  presented  in  ;^>pendix  G. 
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3 .  Antenna  Footprint 

In  order  to  clearly  demonstrate  antenna  coverage,  the 
relative  gain  for  each  beam  was  computed  over  a  5  degree  by  5 
degree  grid  centered  about  the  satellite's  nadir.  The  resvilts 
for  beam  2  are  presented  in  Figure  4.3.  The  beam  pattern  is 
shown  in  terms  of  scan  angle  from  nadir  in  the  in- track  and 
cross-track  directions.  Simileuc  results  were  obtained  for  tha 
remaining  beams  and  are  contained  in  Appendix  G. 
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In  Track 

Aft  Forward 


Relative  dB  Contours 
(Beam  scanned  27  degrees  from  nadir) 


Figure  4.3.  Beaon  2  Relative  dB  Gain  Footprint 
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V.  ALTERNATIVE  DESIGN 


A.  ANIBINA  CJGMFIGDRATION 

An  alternative  design  to  that  presented  in  the  previous 
chapters  is  one  which  tcdces  advantage  of  the  fact  that  the  L- 
band  communications  signal  uses  time  division  multiplexing 
(TM^)  .  The  use  of  TI»1,  as  irrplemented  in  the  payload  for  this 
spacecraft,  means  that  each  of  the  six  antenna  beams  has  a  60 
ms  duty  cycle.  Since  adjacent  pairs  of  antennas  are  not 
operating  at  the  same  time,  the  antenna  could  be  constructed 
such  that  there  are  two  beams  per  antenna  panel  thus  reducing 
the  number  of  panels  required. 

An  integrated  multiple  beam  microstrip  array  similar  to 
that  described  in  [Ref.  27]  cotild  be  constructed.  Two  feed 
lines  would  be  constructed  on  a  lower  svibstrate  which  wo\ild  be 
electromagnetically  coupled  to  the  patch  array  on  the 
sxibstrate  above.  Each  feed  line  would  be  energized  at  the 
appropriate  time  corresponding  to  the  on-time  for  the 
particuleu:  beam.  This  configuration  would  simplify  antenna 
construction  by  eliminating  the  need  for  through-hole  plating 
and  connecting  pins  between  the  feed  network  and  the  patch 
elements  [Ref.  28]  .  The  deployment  mechanism  for  this 
configuration  would  also  be  much  simpler  as  only  two  hinged 
panels  would  have  to  lanfold  as  opposed  to  four.  There  would 
also  be  a  reduction  in  weight  while  performance  would  actually 
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be  enhanced  due  to  reduced  losses  in  the  smaller  feed  network 
required. 

Combining  two  antenna  panels  into  one  and  using  that 
single  panel  for  two  beams  would  almost  double  the  number  of 
patch  elements  available  to  each  beam.  In  addition  to 
increasing  gain,  this  would  allow  better  resolution  in  both 
beam  shape  and  orientation.  Increasing  the  nimnber  of  array 
elements  would  also  eliminate  iinwanted  nulls  which  adversely 
affect  scan  angle. 

B.  MSTTBUtOL  QRIENXATZON 

The  antenna  in  this  study  is  oriented  such  that  individual 
panels  are  perpendicular  to  the  spacecraft '  s  velocity  vector . 
This  orientation  was  chosen  for  stinictural  and  thermal  control 
reasons  during  the  early  stages  of  the  preliminary  design  and 
is  not  optimum.  If  the  design  described  above  were  adopted, 
the  antenna  panels  could  be  oriented  parallel  to  the  velocity 
vector  allowing  cross  track  beam  steering  to  be  accomplished 
by  tilting  two  outside  panels  towards  the  center.  All 
available  patch  elements  could  then  be  dedicated  to  shaping 
the  beam  which  would  result  in  reduced  beam  overlap. 
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VI.  SUNMABY 


A.  ANALYSIS  AND  DESIGN 

Two  different  substrate  materials  were  considered  for  the 
design  of  a  multiple  beam,  phased  array,  microstrip  antenna. 
Design  curves  were  presented  to  illustrate  tradeoffs  between 
weight  and  performance  in  the  substrate  material .  This 
information  was  used  tC'  choose  an  appropriate  substrate 
material  which  would  provide  adequate  performance  at  a  minimum 
cost  in  mass. 

Closed  form  formulae  were  used  to  arrive  at  a  preliminary 
width  and  length  for  microstrip  patch  elements.  A  computer 
model  was  then  used  to  optimize  the  width  and  length 
dimensions .  Array  theory  was  applied  to  derive  patch 
excitation  coefficients. 

B.  PERFORMANCE  EVALOATION 

Patch  excitation  coefficients  were  used  as  input  to  a 
computer  model,  which  employs  modal  expansion  theory,  and 
antenna  radiation  patterns  were  generated. 

While  the  Method  of  Moments  wo\ild  have  yielded  more 
accurate  results,  it  would  have  been  at  the  expense  of 
increased  computer  time.  The  Modal  B:q>ansion  method  used 
provided  adequate  results  with  a  reasonable  amount  of  effort. 
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The  results  demonstrate  that  the  proposed  design  would  not 
meet  all  the  design  goals.  VAiile  the  desired  beam  shapes  were 
obtained,  the  groimd  footprint  of  the  outside  beams  (beam  1 
and  6)  show  that  most  of  the  energy  would  be  concentrated 
outside  the  satellite's  field  of  view. 

C.  C3GNCLDSI0NS 

Analysis  showed  that  this  antenna  could  not  be  constructed 
within  the  allocated  mass  budget.  In  fact,  the  mass  of  the 
sxibstrate  material  alone  would  be  over  two  times  the  allocated 
mass  budget. 

The  use  of  rectangular  microstrip  patches  is  not  an 
efficient  use  of  the  available  array  area.  An  octagonally 
shaped  patch  would  allow  more  elements  to  be  placed  on  a 
panel.  This  would  allow  tighter  control  over  beam  shape  and 
orientation. 

The  alternative  design  discussed  in  C!hapter  IV  would 
address  the  short -comings  in  the  original  design  and  allow  the 
groxrnd  footprint  to  be  entirely  within  the  satellite's  field 
of  view.  However,  operating  two  beams  per  antenna  panel  could 
require  active  phase  shifters,  increasing  cost  as  well  as 
reliability  risk. 

The  use  of  fiber  optic  cables  instead  of  coaxial  transmis¬ 
sion  line  and  etched  copper  feed  networks  would  reduce  both 
weight  and  power  loss. 
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since  the  methods  for  determining  patch  dimensions 
discussed  in  the  literature  all  give  slightly  different 
results,  this  antenna  would  have  to  be  tested  and  trimmed 
after  construction  in  order  to  ensure  the  operating  parameters 
were  met. 

D.  FOTDSE  ^ESIS  OPPOKflUNlTlES 

A  scale  model  could  be  built  and  tested  by  increasing  the 
operating  frequency  to  yield  a  smaller  wavelength.  The 
coupled  feed  model  discussed  in  the  alternative  design  could 
also  be  built  and  tested.  Various  orientations,  beam 
patterns,  and  scaui  angle  could  be  explored  yielding  useful 
data  for  future  phased  array  design  projects. 
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APPBIDIX  A 


LINK  BUDGigr  SUMMARY 


TABLE  A-1.  LINK  EODGEr  [REF.  1] 


1 

User  to 
Satellite 

Satellite 
to  User 

B  Frequency 

1625.0 

NHz 

Frequency 

1625.0 

NHz 

RF  Power 

6 

Watts 

RF  Power 

19 

Watts 

7.8 

dBw 

12.8 

dBw 

Power  Loss 

-1.0 

dB 

Power  Loss 

-2.2 

db 

Antenna  Gain 

3.0 

OBi 

S/C  Ant.  Gain 
(isoflux) 

dBi 

EIRP 

9.8 

dBw 

EIRP 

14.6 

dBw 

Satellite 
Altitude- 
750  NN 

n 

km 

Satellite 

Altitude 

1389.0 

km 

Elevation  Angle 

90.0 

degrees 

Elevation 

Angle 

90.0 

degrees 

Range 

1389.0 

km 

Range 

1389.0 

km 

Free  Space  Loss 

-159.5 

dB 

Free  Space 

Loss 

-159.5 

dB 

-K9.8 

RX  Signal 
Strength 

-155.9 

1  Polarization 

1  Loss 

-0.5 

dB 

Polarization 

Loss 

dB 

Tracking  Loss 

0.0 

dB 

Tracking  Loss 

0.0 

dB 

S/C  Ant.  Gain 
(isoflux) 

A.0-6.0 

dBi 

Antenna  Gain 

mm 

dBi 

RX  Line  Losses 

-1.0 

(S 

RX  Line  Loss 

-0.5 

dB 

L-Band  Aaplifier 
Eff. 

35 

X 

BBmmm 

35 

X 

3.5 

dB 

bjjhshsb 

3.5 

dB 

1  Total  BU 

16.5 

— 

Total  BW 

16.5 

NHz 
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APPENDIX  B 


I.  HASS  ESTIMAIES 

A.  Kevleur  Subtrate 

Density  =  0.05  Ib/in^  Density  =  1.4  gm/cm? 

The  total  siibstrate  voliame  is  given  by: 

1  =  1.4  m  w  =  0.667  m  h  =  0.003  m 
Volume  =lxwxhx2x6 
Volume  =  3.4  X  10*  cm? 

The  total  substrate  mass  is  given  by: 

Mass^gyj,^  =  Volume  x  Density 

^®®iceviar  =  46.5  Kg  Mass^^j,^  =  102.6  lb 

B.  Duroid  Substrate 

Specif ic_gravity  =  2.2  Density  =  2.2  gm/cm* 

Total  sijbstrate  mass: 

MasSoyrojd  “  Volume  x  Density 

=  74  Kg  MasSo„„,j  =  163  lb 
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C.  Copper 

1.  Micjrostrlp  Patches 

Densityj;y  =  8.96  gm/cm? 

Total  nijmber  of  elements  =  9x4x6  =  216 
Element  width  =  5.9  cm  Element  length  =  4.6  cm 

ThicknesSj.y  =  0.0127  cm 

Total  volxjme  of  array  elements  =  length  x  width  x  thicknessj^ 
X  tot  n\]m  elements 

Total  volun'iC  =  74.5  cm* 

Total  Mass,., _ =  Density  x  volume 

Total  Mass.1^^3  =  0.7 

2.  Microstrip  Feed  THnAg 

length^33^  =  1.3  m  width^^^^  =  0.006  m 

Volume^33j  =  length,33j  ^  ^  thickness^y 

Mass^33j  “  Densityc^,  x  Volume^.,^ 

MasSf33j  =  8.9  X  10*^  Kg 

3.  Total  Copper  Mass 

The  total  copper  mass  is  given  by: 

Masses,  =  Mass,i3,^^,  +  Mass^„j 
MasSpy  =  0.701  Kg 

D.  Alianinuni  Ground  Planes 

Since  a  shielded  triplate  feed  is  used,  two  aluminum 
ground  planes  are  needed. 

Density^j  =  2.7  gm/cm*  thickness^^  =  1.0  x  10'^ 

Volume^^  =  1  X  w  X  thickness^j  x  4 
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Mass^i  =  Volume^j  x  Density^^  x  2 
Mass^i  =  1.35  X  10*^  Kg 

II.  SCWkEY 

The  Duroid  substrate  is  approximately  1.5  times  heavier 
than  the  Kevlar  for  a  given  volume.  Therefore,  the  Kevlar 
material  should  be  used  for  the  substrate.  The  total  mass  for 
the  antenna  is  estimated  as  follows: 

Mass  =  Mass^^j,^  +  Massey  +  Mass^i 

A  mass  margin  of  0.01  is  chosen  and  the  total  antenna  mass 
becomes 

^®®antenn«  ®  Mass  X  (1  +  MasS^pjj„) 

=  52  wass„„  .  114.5  lb 
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APPQIDZX  C 


I.  PATCH  DESIGN 

Using  the  transmission  line  model,  the  patch  pcirameters 
are  calculated  as  follows: 

Given:  f  =  1618.25  MHz  =  3.9  h=0.3  cm 

1  =  3(10)Vf  k  =  27cA 

X=18.5cm  k=33.9 


e,+l 


(cm) 


W  =  5.9  cm 


The  effective  dielectric  constant  is  given  by: 

*  2  2  \  W  ) 


=  3.6 


The  patch  length  is  then  calculated  as  follows: 

(e^+0.3)  (-^+0.264) 

AJ  =  0.412h  - -  A1  =  0.1  cm 

(e--0.258)  (-=+0.8) 

®  h 

L  =  - - —  -  2AJ  (cm) 

L  =  4.6  cm 
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II.  GaON 

In  order  to  calculate  gain,  the  directivity  must  be 
calculated  according  to  the  following: 


P  _  _ _ 

|sin  (  )  2tan  (6)  ^sin  (0)  d0 

Losses,  expressed  as  resistance,  must  also  be  foxond  according 

to  the  following  relations: 

Quality  factor:  „  _ _ 

3(10)°ye; 

^  A  f  h 


Radiation  resistance: 
Copper  loss: 

Rc  = 


=  120A,/W 


2.7  (10)-*, 


( 


— )  {-) 
103  w' 


0^ 


Dielectric  loss: 


Rd  = 


81  (10) -3  h  k  Q-i 


The  total  resistance  is  then  Rj  =  R^  +  Rj.  +  Rj  R^  =  379.8 
With  a  of  VSWR  =  2,  Bandwidth  (BW)  and  efficiency  {t\)  are 
found  as  follows: 


BW 


iVSWR-1)  ^^00 
iO  i/VSM) 


T1 


Rj> 


Directivity  and  Gain  are  con:?>uted  as: 

D  =  10  X  log(DQ)  D  =  5.06  dB 

G  =  10  X  log(Tl  X  Dq)  G  =  5.0  dB 
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2 


3 


0.013 

A52.36’ 

0.046 

^-127.6* 

0.046 

A.- 127. 6’ 

0.019 

ii.-57.73" 

0.066 
^12 9. 3* 

0.066 

A.129.3* 

0.024 

Z.-153.8* 

0.083 

L28.18* 

0.083 

iL28.18’ 

0.027 

^103.1" 

0.094 

L-76.9' 

0.094 

A.-76.9’ 

0.028 

iLO.O* 

0.098 

A.180.0* 

0.098 

A.-180.0’ 

0.027 

A.-103.l‘ 

0.094 

A.76.9* 

0.094 
A.76 . 9’ 

0.024 
L153 . 8’ 

0.083 

A-28.18* 

0.083 

A.-28.18* 

0.019 

LSI . 73* 

0.066 

A.- 129. 3’ 

0.066 

A.-129.3’ 

0.013 

ii.-52.36’ 

0.046 

A.127.6’ 

0.046 

A.127.6’ 

0.013 

A.52.36‘ 

0.019 

^-57.73* 

0.024 

^-153.8* 

0.027 

A.103.1* 

0.028 

A.0.0* 

0.027 

A.-103.1* 

0.024 
^153 . 8* 

0.019 

A57.73' 

0.013 

A.-52.36' 
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Element 

Ntunber 


2 


3 


4 


1  ^ 

0.003 
^169. 3* 

0.010 

A.-10.68* 

0.010 

A.-10.68‘ 

0.003 

4.169.3' 

1  ^ 

0.011 

A-55.0‘ 

0.039 

^127.0* 

0.039 

^127.0* 

0.011 
4.-55.  o' 

1  ^ 

0.029 

A.-95.3’ 

0.099 

A.84.6' 

0.099 

A84.6’ 

0.029 

4.- 95. 3' 

4 

0.043 

^-137.7‘ 

0.148 

A.42.3' 

0.148 
^42. 3* 

0.043  1 

4.-137.7'  1 

5 

0.048 

A-180.0’ 

0.167 

/LO.O* 

0.167 

A-O.O’ 

0.048  1 

4.180.0'  1 

6 

0.043 

A.137.7’ 

0.148 

^-42.3’ 

0.148 
^-42. 3* 

0.043  1 

4.137.7'  1 

7 

0.029 

L95 . 3* 

0.081 

A.-84.6’ 

0.099 

A.-84.6' 

0.029  1 

4.95.3'  1 

8 

0.011 
^53 . 0* 

0.039 

A.-127.0’ 

0.039 

A.-127.0* 

0.011  1 
4.53 .  o'  1 

9 

0.003 

/L-169.3’ 

0.010 

A10.68’ 

0.010 

^10.68* 

0.003  1 

4.-169.3'  1 

Element 

Nximber 


0.013 

Ll.l‘ 


0.046 

A.-172.3‘ 


0.046 

A-172.3* 


0.013 

Ll.l' 


0.002 

A.5.8‘ 


0.007 
A- 174. 2* 


0.007 

A.-174.2‘ 


0.002 

15.8* 


0.024 

^-176.1* 


0.083 
A3. 8* 


0.083 
A3. 8* 


0.024 

A-176.1* 


0.051 

A-178.l‘ 


0.176 

A1.9* 


0.176 

A1.9* 


0.051 
A- 178.1* 


0.062 
A- 180.0* 


0.215 

AO.O’ 


0.215 

AO.O* 


0.062 

A180.0* 


0.051 
A178 . 1* 


0.176 

A-1.9‘ 


0.176 

A-1.9’ 


0.051 
A178 . 1* 


0.024 
A176 . 1* 


0.083 

A-3.8‘ 


0.083 

A-3.8* 


0.024 
A176 . 1* 


8 


0.002 

A-5.8* 


0.007 

A174.2* 


0.007 

A174.2* 


0.002 

A-5.8* 


0.013 

A-7.7* 


0.046 

A172.3* 


0.046 

A172.3* 


0.013 

A-7.7* 
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1  Element 
Number 

1 

2 

3 

4 

1 

■siS 

0.046 

^172.3‘ 

0.046 

A172.3’ 

2 

0.002 

^-5.8* 

0.007 

A174.2* 

0.007 

A174.2’ 

0.002 

A-5.8’ 

3 

0.024 
L116 . 1* 

0.083 

^-3.8* 

0.083 

A-3.8’ 

0.024 

A176.1’ 

4 

0.051 

^178.1* 

KBSsM 

HSRSjH 

KEEfll 

0.051 

A178.1’ 

5 

0.062 

^-180.0* 

0.215 

AO.  o’ 

0.215 

AO.  o’ 

0.062 
A180 .  o’ 

6 

0.051 

^-178.1* 

0.176 

A1.9* 

0.176 

A1.9’ 

0.051 

A- 178.1’ 

7 

0.024 

A-176.1’ 

0.083 

A3. 8* 

0.083 

A3. 8’ 

0.024 

A-176.1’ 

8 

0.002 

A5.8’ 

0.007 

A-174.2’ 

0.007 

A-174.2’ 

0.002 

A5.8’ 

9 

0.013 

LI.  7' 

0.046 

A- 172. 3* 

0.046 

A-172.3’ 

0.013 

A7.7’ 
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5.  BEAM  5 


Element 

Nijmber 


0.003 

A.-169.3‘ 


0.011 
LBS .  0* 


0.029 

Z.95.3' 


0.043 

^137.7* 


0.048 

^-180.  o' 


0.043 

A-137.7' 


0.029 

A.-95.3’ 


0.011 
A.- 53.0* 


0.003 

A.169.3’ 


2 

3 

0.010 
^10. 68* 

0.010 

A.10.68’ 

0.039 

A- 127.0* 

0.039 

^-127.0* 

0.099 

A- 84. 6* 

0.099 

A.- 84. 6* 

0.148 

/l-42.3* 

0.148 

A.-42.3* 

0.167 

A.0.0* 

0.167 

A.0.0* 

0.148 

A.42.3* 

0.148 

A.42 . 3* 

0.081 

A84.6* 

0.099 
A.84 . 6* 

0.039 

A127.0' 

0.039 
A.127 . 0* 

0.010 

A-10.68* 

0.010 

A.-10.68* 

0.003 

Z.-169.3' 


0.011 
^55 . 0* 


0.029 

A.95.3* 


0.043 
^13  7. 7’ 


0.048 

A180.0* 


0.043 
^-137. 7* 


0.029 

^-95.3’ 


0.011 

^-53.0* 


0.003 

A.169.3' 


3 


2 


3 


4 


Element 

Number 


0.013 

A.-52.36* 


0.019 

A57.73’ 


0.024 
A.153 . 8* 


0.027 

A-103.1* 


0.028 

6.0.0* 


0.027 
^103 .  l’ 


0.024 
^-153. 8‘ 


0.019 

A.-57.73' 


0.013 

A.52.36‘ 


0.046 

A127.6’ 

0.046 

Z.127.6’ 

0.013 

^-52.36’ 

0.066 

^-129.3* 

0.066 
^-129. 3’ 

0.019 
LSI. 12* 

0.083 
^-28, 18’ 

0.083 

A.-28.18’ 

0.024 
A.153 . 8’ 

0.094 
^76. 9’ 

0.094 

L76 . 9’ 

0.027 

^-103.1’ 

0.098 

A.180 .  o’ 

0.098 

A-180.0’ 

0.028 

A.0.0’ 

0.094 

A.-76.9’ 

0.094 

L-IS.S* 

0.027 
^103 .  l’ 

0.083 

A28.18’ 

0.083 

A28.18’ 

0.024 

^-153.8’ 

0.066 

A.129.3’ 

0.066 
^129. 3’ 

0.019 
^-57. 73’ 

0.046 

A-127.6’ 

0.046 

^-127.6’ 

0.013 

^2.36’ 

APPENDIX  E 


Figure  E>1.  Beam  1  Array  Factor  (^■90') 


Figure  E-2.  Beam  1  Array  Factor  (^«0*) 


Figure  E-8.  Beam  4  Array  Factor  (^sO*) 
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Figure  E-10.  Beam  5  Array  Ftetor  (4>«0*) 
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Figure  E-11.  Beam  6  Array  Factor  (<|>b90*) 


Figure  E-12.  Beam  6  Array  Factor  (4>-0*} 
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APPQ3DIX  F 


I.  DESIGN.FQR 


C.  .  .MIC310STRIP  ANTENNA  DESIGN  CHARACTERISTICS  CALCULATED  BY 
MODAL 

C.  .  .THEORY,  ALL  INPUT  LENGTH  ARE  IN  CENTIMETER,  BY  DR.  HUANG 
COMPLEX  CJ,YWY,AY,DEL(5) ,C1,C2,D4,SKY,FC,WC 
REAL  A,B,EPS 
CJ=(0.,1.) 

PI=3. 14159265 
TPI=2 . *PI 
DPR=180./PI 

TOL= . 0001 

WRITE (6, 11) 

11  FVDRMATC  CENTER  FREQ=?  •  ) 

READ (5, 99)  FREQ 
WRITE (6  1) 

1  FORMAT ('  WIDTH  OF  PATCH=?') 

READ(5,99)A 

99  FORMAT  (F5. 8) 

WRITE (6, 2) 

2  FORMAT ('  RESONANT  LENGTH  OF  PATCH=? ' ) 

READ(5,99)B 

WRIIT!  (6  3 ) 

3  FORMAT ('  DIELECTRIC  CONSTANT*?') 

READ (5, 99) EPS 

WRITE (6, 4) 

4  FORMAT ('  DIELECTRIC  THICKNESS*?') 

READ{5,99)T 

WRITE (6, 5) 

5  FORMAT ('  PROBE  LOCATION*?') 

READ (5, 99) YO 

WRITE (6  6 ) 

6  FORMAT ('  LOSS  TANGENT*?') 

READ(5,99)TLOS 

4  0  FF*1 .  /SQRT  ( 1 .  +10 .  *T/A) 

EPE* (EPS+1 . ) /2 . + (EPS-1 . ) /2 . *FF 

FF*  (A/T+0 .262) / (A/T+0 . 813) * (EPE+0 .3) / (EPE-0 .258) *0.412 
WAMDA*2 . 1*B*SQRT (EPS) 

GW*0 . 00836*A/WAMDA 

BWaO .  01668*FP*A*EPE/WAMDA 

FY*0 . 7747+0 . 5977* (A/B-1 . ) -0 . 1638* (A/B-1 . ) **2 . 
YWY=CMPLX(GW,BW) 
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AY=CJ*TPI*376 . 7/WAMDA*T/A*YWY*FY*B 
DEL(1)  =  (0.,0.) 

DO  20  1=1,4 

C1=2.*AY*(CMPLX(PI,0.)  -DEL(I) ) 

C2=AY*AY+2 . *DEL ( I ) *PI -DEL { I ) *DEL ( I ) -CMPLX ( PI *PI , 0 . ) 
DEL (I+l) =C1/C2- (DEL (I) ) **3 . /3 . 

20  CONTINUE 
D4=DEL(5) 

SKY= (CMPLX (PI , 0 . ) -D4) /B 
FC=CMPLX ( 1 . , 0 . ) -CJ*TLOS 
WC=3 . ElO/SQRT (EPS*FC) *SKY 
FR=REAL (WC) /TPI 
QR=REAL  (WC)  /2  .  /AIMAG  (WC) 

FF=COS (PI*YO/B) 

CC=EPE*8 . 854*1 .E-14*A*B/2 . /T/FF/FF 
RAD=QR/REAL (WC) /CC 

RCU=0 . 00027*SQRT (FR/l .E+9) *B/A*QR*QR 

RDI=30 .  *TLOS*T*WAMDA/EPS/A/B*QR*QR 

RIN=RAD+RCU+RDI 

Q=RIN*REAL  (WC)  *CC 

DFR=FR/Q 

EFF=RAD/RIN 

FR=FR/l.E+9 

DFREQ=FREQ-  FR 

IF(  (ABS(DFREQ) )  .LE.TOL)GO  TO  42 
WRITE ( 6 , 4 1 ) FR , DFREQ 
IP (PR. LT. FREQ)  THEN 
B=B-.0001 
ELSE 

B=B+ . 0001 
END  IF 
GOTO  40 


41  FORMATC  FREQUENCY=' ,F10.5,3X, 'DFREQ=' ,F10.5) 


42  WRITE (6, 31) FR 

31  FORMAT  ( '  RESONANT  FREQUENCY  IN  GIGAHERTZ=  ' ,  FIO  .  5) 
WRITE(6,32)RIN 

32  FORMAT  ( '  INPUT  RESISTANCE= '  ,  FIO .  5 ) 
DFR=DFR/l.E+6/2. 

WRITE  (6  33)  DBT?, 

33  FORMAT  ('  BANDWIDTH  IN  MEGAHERTZ= FIO  .  5) 

WRITE (6  34)EFF 

34  FORMAT ( '  RADIATION  EFFICIENCY= ' , FIO . 5 ) 

PERCBW=1/ (Q*SQRT (2.0)) 

WRITE(6,37)PERCBW 
FORMATC  %  BW=',F10.5) 


37 
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WRITE  ( 6 , 3  8 )  ROJ,  RDI ,  RAD,  Q 

38  FORMAT('RCU=' ,F10.5,2X, •RDI= • , FIO . 5, 2X, 'RAD='  FlO.5, 
2X, 'Q=' ,F10.5) 

WRITE  (6, 39)  A, B 

39  FORMAT ( 'WIDTH=' ,F10. 5, 5X, ' LENGTH= ' , FIO . 5 ) 

STOP 

END 


II.  mCARY.FOR 


C. . .LINEARLY  POLARIZED  PLANAR  ARRAY  FORMED  BY  LINEARLY 
C. . .POLARIZED  MICROSTRIP  ELEMENTS 

C. . .WITH  CROSS-POL,  MODAL  EXPANSION  MULTIMODES  ARE  USED. 
C.  .  .NO  GTD  DIFFRACTION,  ALL  DIMENSIONS  IN  WAVELENGTHS 
C. . .PROGRAMMED  BY  DR.  J.  HUANG  AT  JET  PROPULSION  LAB. 

DIMENSION  DBR(362) ,DBL(182) ,XAX(362) ,XX(2) ,YY(2) , 
1IANG(40) 

DIMENSION  PHASR(182) ,PHASL(182} ,XE(40) ,  YE(40} , 
1AMP(40),  PHAS(40) 

COMPLEX  CJ,CJJ,EX,EY,EZ,XO,YO,ZO,ETH,EPH,ER,EL 
COMPLEX  FX,FY,FZ,X1,Y1,Z1,X2,Y2,Z2,X3,Y3,Z3 
COMMON/DDl/EPS , T, A, B , 12 , IP, XS , YS , ZS , CJJ 
COMMCttJ/DD2 /CJ , PI , TPI , DPR 
COMMON/DD4/THPR, PHPR, RR 

C.  .  .MODIFIED  TO  READ  DATA  FROM  A  GEOMETRY  FILE 
C. . .BY  BILL  BARFIELD,  NPGS 
OPEN (33,  FILE='DATA') 

OPEN (35,  FILE=' GEOMETRY') 

CJ=(0.,1.) 

CJJ=(1.E-15,1.E-15) 

PI=3. 14159265 
TPI*2 . *PI 
DPR=180 . /PI 
WRITE (6  1 ) 

1  FORMAT  ( '  DIELECTRIC  CONSTANT  AND  THICKNESS=?  ' ) 

C  READ (5, 99) EPS, T 

READ (33,*) EPS, T 
WRITE (6, 99)  EPS,T 

99  FORMAT(  F8.6,P8.6) 

98  FORMAT  (18) 

WRITE (6  2 ) 

2  FORMAT  ( '  MICROSTRIP  WIDTH  AND  LENGTH^?  ' ) 
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C  READ(5,99)A,B 

READ(33,*)A,B 
WRITE (6, 99)  A, B 

C. . .A,B  ARE  THE  PHYSICAL  WIDTHS  OF  THE  METALLIC  PATCH 

C . A  ALONE  X-AXIS,  B  ALONE  Y-AXIS,  X-AXIS  IS  HORIZONTAL 

WRITE (6, 3) 

3  FORMAT ( '  FIELD  POINT  DISTANCE  AND  PATTERN  CUT 
/ANGLE=?') 

C  READ (5, 99) RF, PHD 

READ(33,*)RF,PHD 
WRITE  (6, 99) RF,  PHD 
WRITE(6,*)RF,PHD 

WRITE(6,5) 

5  FORMAT ('  FEED  PROBE  DIAMETER  &  DISTANCE  FROM  EDGE=? ' ) 

C  READ(5,99)D,YP 

READ(33,*)D,YP 
WRITE(6,99)D,YP 

WRITE(6,4) 

4  FORMAT  ( '  NUMBER  OF  ELEMENTS^?  ' ) 

C  READ(5,98)N 

READ(35,*)N 
WRITE(6,*)N 

WRITE (6, 6) I 

6  FORMAT ( '  NUMBER ' , 12 , 2X, ' X, Y  LOCATION, AMPLITUDE , PHASE  & 
/ORIENTATION*? ' ) 

DO  105  1=1, N 

C  READ(5,*)XW,YW,AMG,DEG,IANG(I) 

READ  (35 ,  * )  XW,  YW,  AMG, DEG,  lANG  { I ) 

WRITE  ( 6 ,  * )  XW,  YW ,  AMG,  DEG,  lANG  ( I ) 

XE(I)=XW 

YE(I)=YW 

AMP(I)=AMG 

PHAS (I) =DEG/360 . *TPI 

105  CONTINUE 
D=5 . *D 

SK2=EPS*TPI*TPI 
PIB=(PI/B)**2 
PIA=(PI/A)**2 
Fl=l .  /A/B/SK2 
Fl=l  E“7 
CYP=COS (PI*YP/B) 
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F2=2  .  /A/B/  (SK2-PIB)  *CTP 
SD=SIN(PI*D/A)  /  (PI*D/A) 

F3=4 . /A/B/ (SK2-4 . *PIA-PIB) *CyP*SD 
CYP=C»S {TPI*YP/B) 

F4=4 . /A/B/ {SK2-4 . *PIA-4 . *PIB) *CYP*SD 

F5=2  .  /A/B/  (SK2-4  .  *PIB)  *CYP 

F6=2  .  /A/B/  (SK2-4  .  *PIA)  *SD 

F7=2 .  /A/B/ (SK2-  9 . *PIB) *COS (3 . *PI*YP/B) 

C  WRITE(6,7)F1,F2,F3,F4 

7  FORMAT(  'Fl='  ,F10.5,5X,  ' F2= ' , FIO . 5, 5X,  ' F3= ' , FIO  . 5, 5X,  ' 
/F4=' ,F10.5) 

C  WRITE(6,8)F5,F6,P7 

8  FORMATC  F5=',F10.5,5X, •F6=',F10.5,5X, 'F7=',F10.5) 

DO  200  1=1,181 

THDD=FLQAT ( I - 1 ) /2 . -45 - 
THD=ABS  (THDD) 

THR=THD/DPR 

PHR=PHD/DPR 

IF (THDD . LT . 0 . ) PHR=PHR+PI 
XS=RF*SIN(THR) *COS (PHR) 

YS=RF*SIN(THR)  *SIN(PHR) 

ZS=RF*COS  (THR) 

FX=CJJ 

FY=CJJ 

FZ=CJJ 

DO  104  14=1, N 
EX=CJJ 
Ey=CJJ 
EZ=CJJ 

DO  102  12=1,4 
XA=  (B+T) /2 . 

IF(I2.EQ.2)XA=-XA 
IF(I2.EQ.3)XA=(A+T)  /2 . 

IF(I2.EQ.4)XA=-(A+T)  /2 . 

X=XS-XE(I4)-XA 

Y=yS-YE(I4) 

Z=ZS 

IF (12 . LT . 3 ) X=XS-XE ( 14 ) 

IF(I2.LT.3)Y=YS-YE(I4) -XA 
RR=SQRT (X*X+Y*Y+Z*Z) 

XU=X/RR 

YU=Y/RR 

ZU=Z/RR 

IF(ABS(XU)  .LT.l.E-8)XU=l.E-8 
IF(ABS(YU) .LT.l.E-8)YU=l.E-8 
PHPR=ATAN2  (YU,  XU) 

IP ( PHPR . LT . 0 . ) PHPR=PHPR+TPI 
IF ( 12 . LT . 3 ) PHPR=PHPR-Pl/2 . 

IF ( PHPR . LT . 0 . ) PHPR=PHPR+TPI 
THPR=AODS  (ZU) 
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CALL  SLOT4(XO,YO,ZO,X1,Y1,Z1,X2,Y2,Z2,X3,Y3,Z3) 

IF(IANG(I4) .EQ.270)GO  TO  50 

IF(IANG(I4) .EQ.90)GO  TO  60 

IF(IANG(I4) .EQ.180)GO  TO  70 

IF(I2.EQ.2)GO  TO  41 

IF(I2.EQ.3)GO  TO  42 

IF(I2.EQ.4)GO  TO  43 

EX=EX- {F1-F2+F5-F7) *Y0- (F4-F3+F6) *Y2 
EY=EY+ (F1-F2+F5-F7) *X0+ (F4-F3+F6) *X2 
EZ=EZ+ (F1-F2+F5-F7) *Z0+ (F4-F3+F6) *Z2 
GO  TO  102 

41  EX=EX+ (F2+F1+F5+F7) *Y0+ (F3+F4+F6) *Y2 
EY=EY- (F2+F1+F5+F7) *X0- (F3+F4+F6) *X2 
EZ=EZ- (F2+F1+F5+F7) *Z0- (F3+F4+F6) *Z2 
GO  TO  102 

42  EX=EX+ (F1-F6) *X0+ (F2-F3) *X1+ (F4-F5) *X2-F7*X3 
EY=EY+ (F1-F6) *Y0+ (F2-F3) *Y1+ (F4-F5) *Y2-F7*Y3 
EZ=EZ+ (F1-F6) *Z0+ (F2-F3) *Z1+ {F4-F5) *Z2-F7*Z3 
GO  TO  102 

43  EX=EX- (F1-F6) *X0+ (F3-F2) *X1- (F4-F5) *X2+F7*X3 
EY=EY- (F1-F6) *Y0+ (F3-F2) *Y1- (F4-F5) *Y2+F7*Y3 
EZ=EZ- (F1-F6) *Z0+ (F3-F2) *Z1- (F4-F5) *Z2+F7*Z3 
GO  TO  102 

50  IF(I2.EQ.2)GO  TO  51 
IF(I2.EQ.3)GO  TO  52 
IF(I2.EQ.4)GO  TO  53 

EX=EX- ( (F1-F6) *Y0+ (F3-F2) *Y1+ (F4-F5) *Y2+F7*Y3) 
EY=EY+ ( (F1-F6) *X0+ (F3-F2) *X1+ (F4-F5) *X2+F7*X3) 
EZ=EZ+ ( (F1-F6) *Z0+ (F3-F2) *Z1+ (F4-F5) *Z2+F7*Z3) 
GO  TO  102 

51  EX=EX- (- (F1-F6)*Y0+{F2-F3)*Y1- (F4-F5) *Y2-F7*Y3) 
EY=EY+ (- (F1-F6) *X0+ {F2-F3) *X1- {F4-F5) *X2-F7*X3) 
EZ=EZ+ (- (F1-F6) *Z0+ (F2-P3) *Z1- {F4-F5) *Z2-F7*Z3) 
GO  TO  102 

52  EX=EX+ ( (F1-F2+F5-F7) *X0+ (F4-F3+F6) *X2) 

EY=EY+ ( (F1-F2+F5-F7) *Y0+ (F4-F3+F6) *Y2) 

EZ=EZ+ ( (F1-F2+F5-F7) *Z0+ (F4-F3+F6) *Z2) 

GO  TO  102 

53  EX=EX+( (-F2-Fl-F5-F7)*X0+(-F3-F4-F6)*X2) 

EY=EY+ ( (-F2-F1-F5-F7) *Y0+ (-F3-F4-F6) *Y2) 

EZ=EZ+ ( (-F2-F1-F5-F7) *Z0+ (-F3-F4-F6) *Z2) 

GO  TO  102 

60  IF(I2.EQ,2)GO  TO  61 
IF(I2.EQ.3)GO  TO  62 
IF(I2.EQ.4)GO  TO  63 

EX=EX- ( (F1-F6) *Y0+ (F2-F3) *Y1+ (F4-F5) *Y2-F7*Y3) 
EY=EY+ (F1-F6) *X0+ (F2-F3) *X1+ (F4-F5) *X2-F7*X3 
EZ=EZ+ (F1-F6) *Z0+ (F2-F3) *Z1+ (F4-F5) *Z2-F7*Z3 
GO  TO  102 

61  EX=EX+ (F1-F6) *Y0+ (F2-F3) *Y1+ (F4-F5) *Y2-F7*Y3 
EY=EY- ( (F1-F6) *X0+ (F2-F3) *X1+ (F4-F5) *X2-F7*X3) 
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EZ=EZ- ( {F1-F6) *Z0+ (F2-F3) *Z1+ (F4-F5) *Z2-F7*Z3) 

GO  TO  102 

62  EX=EX+ (F1+F2+F5+F7) *X0+ (F3+F4+F6) *X2 
Ey=EY+ (F1+F2+F5+F7) *Y0+ (F3+F4+F6) *Y2 
EZ=EZ+ (F1+F2+F5+F7) *Z0+ (F3+F4+F6) *Z2 
GO  TO  102 

63  EX=EX+ (F2-F1-F5+F7) *X0+ (F3-F4-F6) *X2 
EY=EY+ (F2-F1-F5+F7) *Y0+ (F3-F4-F6) *Y2 
EZ=EZ+ (F2-F1-F5+F7) *Z0+ (F3-F4-F6) *Z2 
GO  TO  102 

70  IF(I2.EQ.2)GO  TO  71 
IF(I2.EQ.3)GO  TO  72 
IF(I2.EQ.4)GO  TO  73 

EX=EX- ( (F1+F2+F5+F7) *Y0+ {F3+F4+F6) *Y2) 

EY=EY+ (F1+F2+F5+F7) *X0+ (F3+F4+F6) *X2 
EZ=EZ+ (F1+F2+F5+F7) *Z0+ (F3+F4+F6) *Z2 
GO  TO  102 

71  EX=EX- ( (F2-F1-F5+F7) *Y0+ (F3-F4-F6) *Y2) 

EY=EY+ (F2-F1-F5+F7) *X0+ (F3-F4-F6) *X2 
EZ=EZ+ (F2-F1-F5+F7) *Z0+ (F3-F4-F6) *Z2 
GO  TO  102 

72  EX=EX+ (F1-F6) *X0+ {F3-F2) *X1+ (F4-F5) *X2+F7*X3 
EY=EY+ (F1-F6) *Y0+ (F3-F2) *Y1+ (F4-F5) *Y2+F7*Y3 
EZ=EZ+(F1-F6) *Z0+(F3-F2) *Z1+(F4-F5) *Z2+F7*Z3 
GO  TO  102 

73  EX=EX+ (F6-F1) *X0+ (F2-F3) *X1+ (F5-F4) *X2-F7*X3 
EY=EY+ (F6-F1) *Y0+ (F2-F3) *Y1+ {F5-F4) *Y2-F7*Y3 
EZ=EZ+ (F6-F1) *Z0+ (F2-F3) *Z1+ (F5-F4) *Z2>F7*Z3 

102  COKTINUE 

EX=FX+EX*AMP ( 14 ) *CEXP ( -CJ*PHAS ( 14 ) ) 

FY=FY+EY*AMP ( 14 ) *CEXP ( -CJ*PHAS ( 14 ) ) 

FZ=FZ+EZ*AMP { 14 ) *CEXP ( -CJ*PHAS ( 14 ) ) 

104  CONTINUE 

ETH=FX*COS  (PHR)  *COS  (THR)  +PY*SIN  (PHR)  *COS  (THR)  -FZ*SIN 
/(THR) 

EPH=-FX*SIN  (PHR)  +FY*COS  (PHR) 

C  ER= (ETH+CJ*EPH) /SQRT ( 2 . ) 

C  EL= (ETH-CJ*EPH) /SQRT ( 2 . ) 

ERsETH 
EL^EPH 
AER=CABS (ER) 

AELsCABS (EL) 

IF (AER . LT . 1 . E- 15 ) AER=1 . E- 15 
IF (AEL . LT . 1 . E- 15 ) AEL=1 . E- 15 
IF (AER . LT . 1 . E- 8 ) ER=CMPLX ( 1 . , 1 . ) 

IF (AEL . LT . 1 . E-8) EL=CMPLX (1 . , 1 . ) 

DBR (I) >20 . *ALOG10 (AER) 

DBL(I) >20 . *ALOG10 (AEL) 

PHASR  ( I )  >ATAN2  (AINAG  (ER) ,  REAL  (ER) )  *DPR 
PHASL  ( I )  >ATAN2  (AIMAG  (EL)  ,  REAL  (EL) )  *DPR 
XAX(I)=THDD 
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nonnoonn  on 


c 


WRITE (6 , * ) XAX, DBR, DBL 


200  CONTINUE 

YMAXL=-1000. 

YMAXR=-1000. 

DO  22  1=1  181 

IF ( YMAXR . LT . DBR ( I ) ) YMAXR=DBR { I ) 

IF ( YMAXL . LT . DBL ( I ) ) YMAXL=DBL ( I ) 

22  CONTINUE 

YMAX=YMAXR 

IF  (YMAXR .  LT .  YMAXL)  YMAX=YMAXL 
DO  21  1=1,181 
DBR(I)=DBR(I) -YMAX-2. 

DBL(I)=DBL(I) -YMAX-2. 

IF(DBR(I) .LT.-40.)DBR(I)=-40. 

IF (DBL (I) .LT. -40.) DBL (I) =-40. 

21  CONTINUE 

CALL  BGNPLT 

CALL  PLFORMCLINLIN' ,7.43,4.6) 

XX(1)=-180. 

XX(2)=180. 

YY(l)=-40. 

YY(2)=0. 

CALL  PLSCAL (XX, 2, 080908, YY, 2, 040504) 

CALL  PLGRAF( 'LP  MICROSTRIP  ARRAY' , 'ANGLE' , 'DB' ) 
CALL  PLCURV(XAX,DBR,181,0,0) 

CALL  PLNUP 
CALL  PLNTYP(5) 

CALL  PLNDN(0.,0.) 

CALL  PLCURV(XAX,DBL,181,0,0) 

CALL  ENDPLT 

J=181 
T=-181 

DO  211  1=182,362 
XAX(I)=(I-l)/2.-45. 

IF(XAX(I) .GT.90.)XAX(I)=XAX(I)-180. 

DBR(I)=DBR(J) 

J=J-1 
211  CONTINUE 

OPEN(38,  FILE=' RESULTS. PRN') 

WRITE (6, 699) 

699  FORMAT ( '  WRITING  RESULTS  TO  FILE. . . ' ) 

700  FORMAT (F8. 3, IX, FS. 3,  IX) 

DO  701  1=1,361 
WRITE(38,*)XAX(I) ,DBR(I) 


€B 


701 


CONTINUE 


WRITE (6,*)  '  DONE. . .  ' 

1000  CLOSE (33) 

CLOSE (35) 

CLOSE (38) 

STOP 

END 

SUBROUTINE  SLOT4  (EX,EY,EZ,EX1,EY1,EZ1,EX2, 
/EY2 , EZ2 , EX3 , EY3 , EZ3 ) 

C. . .RADIATION  FROM  DIFFERENT  MODE  SLOTS 

COMPLEX  CJ,FAC,EPH,ETH,EX,EY,EZ,CJJ 

COMPLEX  EXl,Eyi,EZl,EX2,EY2,EZ2,EX3,EY3,EZ3 

COMMON/DDl/EPS , T , AA, BB , 12 , IP , XS , YS , ZS , CJJ 

COMMON/DD2/CJ, PI , TPI , DPR 

C0MM0N/DD4/THR,  PHR,  R 

FAC=CEXP ( -CJ*TPI*R) /R 

A=T/2 . 

B= (AA+T) /2 . 

IF ( 12 . GT . 2 ) B= (BB+T) /2 . 

CPH=COS(PHR) 

SPH»SIN(PHR) 

CrH=COS (THR) 

STH=SIN(THR) 

ARG1=TPI*A*CPH*STH*SQRT (EPS) 
ARG2=TPI*B*SPH*STH 
IF(ABS(ARG1) .LT.l.E-4)GO  TO  11 
F1=SIN(ARG1)  /ARGl 
GO  TO  12 

11  Fl=l . 

12  IF(ABS(ARG2) .LT.l.E-4)GO  TO  13 
F2=SIN (ARG2 ) /ARG2 

GO  TO  14 

13  F2=l . 

14  IF(ABS(ABS(ARG2)-PI/2.) .LT.l.E-4)GO  TO  15 
F3=ARG2*COS  (ARG2)  /  (ARG2*ARG2-PI*PI/4  . ) 

GO  TO  16 

15  F3=-1./PI*ARG2 

16  F4=COS (ARG2 ) / (ARG2  *ARG2 -PI*PI ) 
IF(ABS(ABS(ARG2)-1.5*PI) .LT.l.E-4)GO  TO  17 
F5=ARG2*COS (ARG2) / (ARG2*ARG2-PI*PI*9 . /4 . ) 

GO  TO  18 

17  F5=-l./2. 

18  EPH=-CJ*F1*F2*SPH*CTH*FAC 
ETH=CJ*P1*F2 *CPH*FAC 
EX=ETH*CTH*CPH-EPH*SPH 
EY=ETH*CTH*SPH+EPH*CPH 
EZ=-ETH*STH 
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EPH= - FI *F3  *SPH*CTH*FAC 
ETH=F1*F3  *CPH*FAC 
EX1=ETH*CTH*CPH-EPH*SPH 
EY1=ETH*CTH*SPH+EPH*CPH 
EZ1=-ETH*STH 

EPH=CJ*PI*F1*F4*SPH*CTH*FAC 

ETH=-CJ*PI*F1*F4*CPH*PAC 

EX2=ETH*CTH*CPH-EPH*SPH 

EY2=ETH*CTH*SPH+EPH*CPH 

EZ2=-ETH*STH 

EPH=F1*F5*SPH*CTH*FAC 

ETH= - F1*F5 *CPH*FAC 

ex3=eth*c:th*cph-eph*sph 

EY3 =ETH*CTH*SPH+EPH*CPH 
EZ3=-ETH*STH 
RETURN 
END 


III.  GBGM.FQR 

C. . .ROUTINE  TO  CREATE  GEOMETRY  FOR  MICARY.FOR 

DIMENSION  X(4) ,Y(9) ,SMAGD(9,4) ,SPHASED(9,4) 
REAL  L,W,AMP,  PHASE,  LAMBDA,  FREQ,  SPACE 

OPEN (5 , FILE= ' GEOMETRY ' ) 

OPEN(33,FILE='MAG') 

OPEN (35 , FILE= » PHASE ' ) 

C.  .  .READ  IN  MAGNITUDE  AND  PHASE  DATA  FROM  FILE 


NUMX=4 

NUMY=9 

NUMTOT=NUMX*NDMY 


DO  4  I=1,NUMY 

READ(33,*)SMAGD(I,1) ,SMAGD(I,2) ,SMAGD(I,3) ,SMAGD(I,4) 
READ(35,*)SPHASED(I,1) ,SPHASED(I, 2) , SPHASED(I, 3) , 
/SPHASED  (1,4) 

WRITE(6,*)^GD(I,1)  ,SMAGD(I,2)  ,SMAGD(I,3)  ,SMAGD(I,4) 
WRITE (6,*) SPHASED (1,1) , SPHASED ( I , 2 ) , SPHASED ( I , 3 ) , 
/SPHASED  (1,4) 

4  CONTINUE 

12  FORMAT(F10.5,F10.5,F10.5,F10.5) 

C. . .FREQUENCY  IN  HZ 
FREQ=1618.25E6 
LAMBDA=3 .E8/FREQ 
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C. . .DIMENSIONS  IN  CM 
L=4.44 
W=5.9 
AMP=1.0 
PHASE= . 0 
ORIENT=0 . 

SPACE= . 55*LAMBDA 

C. . .CONVERT  TO  WAVELENGTHS 
L=L/100 . 

L=L/LAMBDA 
W=W/100 . 

W=W/LAMBnA 

DO  10  I=1,NUMX 
X (I) = (I-l) *SPACE+ (I-l) *W 
10  CONTINUE 

DO  15  I=1,NUMY 
Y(I)  =  (I-1)*SPACE+(I-1)*L 
15  CONTINUE 

NUMEL=NUMX*NUMY 

WRITE(5,98)NUMEL 

98  FORMAT  (13) 

K=0 

N=0 

ORIENT* 0 
DO  20  1=1,9 
DO  22  J=l,4 

K=K+1 

WRITE(5,99)X(J)  ,Y(I)  ,SMAGD(I,J)  , SPHASED ( I ,  J)  ,N 

22  CONTINUE 
20  CONTINUE 

99  FORMAT  (F8 . 3 ,  IX,  F8 . 3 ,  IX,  F8 . 3 ,  IX,  F8 . 3 ,  IX,  II) 

STOP 

END 
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APPm>IX  6 


I.  RADIATION  PATTERNS 


Figure  G-1.  Beam  1  Radiatioai  pattern 


Figure  Q-2.  Beam  2  Radiatioai  Pattern  ((i>si90*’) 
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Relative  Gain  (dB)  j  |  Relative  Gain  (dB) 


Figure  G-3.  Beam  3  Radiation  Pattern  (i^s90^} 


Figuire  0-4.  Beam  4  Radiation  Pi.ttem  (^*90**) 
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Figure  G>9.  Siinultaneous  Radiation 

III.  ANIENNA  FOOTPRINT 


Figure  0-10.  Beam  1  Relative  dB  Gain  Footprint 
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Figure  6-13.  Beam  4  Relative  dB  Gain.  Foot^nrint 
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Figure  G-14.  Beam  5  Relative  dB  Gain  Footprint 
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Monterey,  CA  93943-5121 

6.  LT  Jcimes  A-  Frank  1 

8508  Blue  Bird  Woods  Ct. 

Lorton,  VA  22079 

7.  LCDR  Randall  Butler  1 

PSC  477  Box  417 

FPO  AP  96306-1299 

8.  LT  Stephanie  Kessler  1 

COMUSNAVCENT  Code  N2 

FPO  AE  09501-6009 

9.  LT  William  L.  Barfield  2 

2836  Kelly  Square 

Vienna,  VA  22181 
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